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We report the first experimental realization of single-qubit manipulation for single spinwaves 
stored in an atomic ensemble quantum memory. In order to have high-fidelity gate operations, we 
make use of stimulated Raman transition and controlled Lamor precession jointly. We characterize 
the gate performances with quantum state tomography and quantum process tomography, both 
of which imply that high-fidelity operations have been achieved. Our work complements the ex¬ 
perimental toolbox of atomic-ensemble quantum memories by adding the capability of single-qubit 
manipulation, thus may have important applications in future scalable quantum networks. 

PACS numbers: 42.50.Dv 


Many physical systems have been experimentally stud¬ 
ied for quantum information applications pQ. Among 
them, atomic ensembles [2] are mainly famous for the ca¬ 
pability of long-term storage of quantum states and the 
capability of efficient interaction with single-photons. In 
recent years, significant achievements on memory lifetime 
El 0 and efficiency El E] have been made. Neverthe¬ 
less, direct manipulation of single qubits in an atomic- 
ensemble quantum memory has not been realized so 
far, which limits further applications of atomic-ensemble 
quantum memories. For instance, in order to teleport a 
qubit from one ensemble to another [5] or entangle re¬ 
mote atomic ensembles |7j through interfering indepen¬ 
dent photons, random single-qubit manipulations on one 
ensemble are required in order to restore the target states 
conditioned on different Bell state measurement results. 
However, due to the incapability of qubit manipulation, 
in all previous experiments 00, one has to convert the 
spinwave states into single-photon states and apply the 
qubit rotations on photons instead. Inefficiency is the 
main drawback for this alternative way. For large-scale 
applications, the overall efficiency can be extremely low. 
If the spinwave states can be directly manipulated and 
detected in the atomic ensemble quantum memories, the 
efficiency will be largely enhanced. 

Arbitrary single-qubit rotation is an essential element 
in quantum information science, and has been success¬ 
fully implemented in single emitter systems, such as sin¬ 
gle neutral atoms ED], single ions m and single quan¬ 
tum dot QH 01 and single NV center mm- However, 
coherently manipulate the spinwave qubit in atomic en¬ 
sembles remains challenging. For single-emitter systems, 
single-qubit manipulation can be easily realized by ad¬ 
dressing the single emitter using microwaves or radio fre¬ 
quency fields. Nevertheless, qubit manipulation for a sin¬ 
gle spinwave qubit requires coherently manipulate all the 
atoms simultaneously, since a single spinwave is a collec¬ 
tive superposition of all the atoms with only one excita¬ 
tion dim. The main difficulty of manipulating such 


a single spinwave lies in how to manipulate a single col¬ 
lective excitation mode precisely and coherently without 
influencing (depopulating or exciting) the rest majority 
atoms. Even a single indeliberate excitation from the ma¬ 
jority of atoms could possibly ruin a stored single spin- 
wave. The manipulation process has also to preserve the 
spinwave wave-vector since its amplitude and direction 
determines the mode direction for the retrieval photons. 
Preliminary study on manipulation of classical spinwaves 
has been performed previously by H. Wang et al. Either 
by manipulating the majority atoms [T5| or the excited 
atoms 01 , they observed the flopping between two spin- 
wave modes. Nevertheless, their experiments still lie in 
the classical regime, and number of excitation is on the 
order of ~ 10 6 which is far more larger than a single 
spinwave. 

In this paper, we report the first experiment of single¬ 
qubit manipulation for single spinwaves in an atomic 
ensemble quantum memory. High-fidelity operation is 
achieved by making use of stimulated Raman transition 
and controlled Larrnor procession jointly. An arbitrary 
single spinwave state is heraldedly prepared through the 
process of spontaneous Raman scattering and making 
projective measurement on the scattered single photons. 
For the realization of an arbitrary single-qubit operation, 
the capabilities of rotating an arbitrary angle along two 
orthogonal axes in the Bloch sphere are required. For 
the rotation along the axis in the direction of poles, we 
make use of controlled Larrnor procession. For the rota¬ 
tion along an axis in the equatorial plane, we make use 
of stimulated Raman transitions. Experimental results 
are characterized with quantum state tomography and 
quantum process tomography. 

Our experimental scheme is shown in Fig. [l] An en¬ 
semble of 87 Rb atoms are captured using a standard 
magneto-optical trap (MOT). With 2 ms polarization 
gradient cooling, the atomic ensemble is cooled to a tem¬ 
perature of about 10 /jK. All the atoms are initially 
pumped to the state of |g) = |F = 1 ,m/ = 0). To 
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create a single spinwave exication |16j , as shown in Fig. 
B a )’ a a polarized write pulse is applied to couple the 
transition of | g) -> |e) with |e) = | F' = 2 ,mp = — 1) 
weakly. Raman scattered a + and a~ idler photons are 
selected through spacial and frequency filtering. A er + 
idler photon heralds the creation of a single excitation 
in | J,) = |F = 2 ,to/ = 2), while a cr~ idler photon 
heralds the creation of a different single excitation in 
| f) = l-F = 2, to/ = 0). For each write pulse, the 
heralding probability is 3 x 10~ 3 , which guarantees that 
probability for creating multiple (> 2) excitations is neg¬ 
ligibly small. Interference between these two channels 
gives rise to the entanglement between an idler photon 
and an atomic spinwave [201 Eli i n the form of 


a Write | Rotation Read 




I’J'ap) = \/2jl\s i )\<T + ) - \/375|s t )|cr ), 

where the amplitude and relative phase are determined 
by the atomic transitions and the collective spin states 
are defined as |s^) = Hje lk, ' Ti \g\... ij ...) and (s/) = 

J^je tk ’' ri \gi... tj with k s = k w — fcj. Once the 

idler photon is measured under specific polarizations, the 
spinwave will be projected to the corresponding states. 

In order to prepare an arbitrary single spinwave state 
of \ip s ) = cos 0 |s 4 .) + sin0e l ^|s-|-), one just needs to 
project the idler photon onto a corresponding state of 
\4>i) = {\/5/5cos9— \/2/5 sin 9e~ l< ^)\H}+i{^/3/5 cos 9 + 
a/ 2/5 sin de - ^)^), where a normalizing factor is omit¬ 
ted. The projection of idler photon can be done through 
standard linear-optics method [22]. Before performing 
qubit rotations on this spinwave state, we first charac- 
tize the state preparation fidelities. The spinwave state 
is measured by applying a strong read pulse which con¬ 
verts |sj,} to a a~ polarized signal photon and converts 
|s-(-) to a + polarized signal photons. Arbitrary single¬ 
qubit measurement is realized through measuring the sig¬ 
nal photon in arbitrary polarization bases, which is also 
one of the highlights of our system. In our experiment, 
we select six spinwave states to prepare: |s/, |s-f), |s_d), 
|sa), | Sr) and | s L ), with \s D/A ) = l/v^dsj.) ± |s T )) 
and | s R / L ) = l/V2{\ Si ) ± *| s t))- We ma ^e use of 
quantum state tomography and the maximum likelihood 
method[23J to characterize the prepared states, and cal¬ 
culate the fidelities in comparison with the target states. 
Results are shown in Tab]!] and an average fidelity of 
97.2% is achieved. Standard errors for each state are cal¬ 
culated by poissonian sampling of the measured counts 
for 200 times. 

In order to realize arbitrary single-qubit manipulation, 
the ability of rotation along two orthogonal axes in the 
Bloch sphere are required. As shown in Fig. [ljb) , we first 
make use of the Larmor procession process to realize the 
rotation of !?-(</>) which corresponds to rotation along an 
axis in the direction of poles. As the magnetic moment 
of the |s/} state is zero, the Larmor phase of this state is 
constant during the experiment. While for the |sj.) state, 
its magnetic moment precesses around the bias magnetic 
field, with an angluar frequency of wp = 2nBo x 1.4 



FIG. 1: Experiment scheme for arbitrary rotations of spin- 
wave qubit, a. Energy levels: The write process projects the 
spinwave onto a superposition state between |sj.) and|s-|-}, by 
measuring the idler photon along a specific polarization state. 
Then a series of rotations between two basis states are carried 
out. At last, the spinwave qubit state is verified by convert¬ 
ing it back to a signal photon during the read process. The 
levels chosen here are designed to provide almost identical 
retrieval efficiencies for both spinwave states, since the tran¬ 
sition strengths of | F' = 2 , m/ = —1/ + 1} —> \g) are the 
same, b. Two elemental rotations in the Bloch sphere. One 
is the rotation along z axis, which is due to different Larmor 
precession frequencies of these two spin states. The other is 
along an axis in the x-y plane, where the orientation angle is 
determined by the relative phase between two Raman beams, 
c, Experimental Layout. The write beam has a waist of 150 
pm, a power of 1.5 pW and a detuning of -10 MHz. The idler 
photon is collected in the direction with an angle of 1.5° rel¬ 
ative to the write beam, and has a mode waist of 100 pm. It 
is further frequency-filtered with a Fabry-Perot cavity. While 
the signal photon is filtered by a Rubidium vapor filter to 
avoid back reflections into the idler channal during write pro¬ 
cess. The polarization of a single Raman beam is prepared 
to provide both helicities for stimulated Raman transitions 
between . 


MHz/Guass. Thus the relative phase between Is^/s-f) 
will evolve as 4>{t) = u>Lt, which effectively is a rotation 
around z-axis in Bloch sphere. A specific rotation angle 
can be realized by changing the duration of the Larmor 
procession. With the bias field applied, we measured the 
Larmor frequency to be wp = 27T x 180kHz. The rotation 
matrix of Larmor operation is, 
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Idler State 

Fidelity(%) 

Spinwave State 

Fidelity(%) 

l*> 

98.3(1.6) 

K) 

97.6(1.0) 

\y) 

97.4(1.2) 

h) 

99.0(0.6) 

\ D ) 

96.2(1.2) 

N 

95.1(1.3) 

k> 

97.6(1.1) 

N 

96.3(1.3) 

\ R ) 

97.4(1.1) 

\Sr) 

97.5(1.3) 

l*> 

99.1(0.5) 

\S L ) 

97.8(1.2) 


TABLE I: Fidelities of state preparations. (Left) Fidelities be¬ 
tween measured and expected signal photon states for differ¬ 
ent idler photon polarizations. Fidelities are relatively higher 
due to the polarization needed is more accurate to prepare. 
(Right) Fidelities between measured and expected signal pho¬ 
ton states for different spinwave states. 


where is the duration of Larmor operation. 

We make use of stimulated Raman transition to re¬ 
alize the rotation along an axis in the equatorial plane. 
The detailed scheme and experimental arrangement are 
shown in Fig. 0 First, we notice that if both Raman 
beams are red (or blue) detuned relative to the \F' = 2) 
and | F' = 1) hyperfine states of the D1 line, destructive 
interference between these two Raman transition chan¬ 
nels will occur. To avoid this, we set the frequency of the 
Raman beams to be at the middle of these two states, 
as shown in Fig. 0a), where the detuning is A « 408 
MHz. Thus these two Raman transitions constructively 
interfere with each other, and lead to faster rotation op¬ 
erations with limited laser power. Second, it’s found that 
when the Raman beams are applied, they couple not only 
the transition of |s-j-) •<+ |sj.) but also an unexpected tran¬ 
sition of |sf) e■» |s aux ) = \F = 2, rriF = +2). As the state 
Isaux) cannot be converted to signal photons in the read 
process, we should minimize the population leakage into 
this state, and restrict the Rabi flopping to occur mainly 
between Is^.) and |s-|-). 

The Raman beam in our setup has a total power 
of about 7 mW, with a waist of 1.9 mm. The po¬ 
larization state of the beam is adjusted to be \i/jr) = 
\] l/7|rr + ) + e l(pR y/Wj7\a~), where ct+/( 7~ terms serve 
as the k + /k~ Raman beams in Fig. 0a) respectively, 
and (pa is the relative phase between these two terms. 
With such a polarization, AC Starks shifts due to the 
Raman beams of |s-|-) and |s aux ) states are calcu¬ 
lated to be about +40 kHz, -140 kHz and +240 kHz, 
respectively. Thus the frequency splitting between | sq-) 
and |sj.) induced by AC Starts effect is -180 kHz, which 
just cancels the Zeeman splitting (+180 kHz) between 
these two states. However, the overall frequency split¬ 
ting between |Sauic) and |s-|-) are shifted to an even larger 
value of 560 kHz. The theoretical two-photon Rabi fre¬ 
quencies of |s-|-) <+ |s|) and |s-f) <+ |s aux .) are both 240 
kHz. As the two-photon detuning of the |s-j-) +»• \s aux ) 
transition is significantly larger than the Raman Rabi 


frequency, the maximum population leakage onto |s aua; ) 
during Rabi flopping is estimated to be ~ 1%. Thus 
we can restrict the population mainly flops between our 
target spinwave states of |s-|-) and |sj,). 

Then we can successfully create a pure rotation along 
a definite axis in the x-y plane of the Bloch sphere when 
a Raman pulse is applied. Note the angle <p between the 
Raman rotation axis and x axis is determined by the rel¬ 
ative phase <pr between k + and k~ Raman components. 
An additional merit of our scheme is that this phase can 
be precisely controlled and it’s stable in itself, excluding 
other electronic phase stabilizing controls. The rotation 
matrix of Raman operation can be written as, 

, \ _ f" cos <j>(t R ) -ie %VR sin (j>(t R ) 

R > —ie~ llpR sin4>(tR) cos ^(t R ) 

where 4>{tR) = f2#t#/2, and tR is the duration of applied 
Raman pulse. By choosing p>R = 0, rotation along x-axis 
can be realized. While with ip R = —7t/ 2, the rotation is 
along y-axis, as shown in Fig. 0b)- With Raman beam 
paramters stated above, the measured Raman Rabi fre¬ 
quency is about 190 kHz. The deviation with theoretical 
predications may be due to the imperfect interference be¬ 
tween two Raman transition channels. 

We first evaluate the gate performances with quan¬ 
tum state tomography [23] by measuring the quantum 
state of the spinwave after a single-qubit operation and 
compare it with the ideal target state using F(pi,p 2 ) = 
{tr[( x /pi/ 02 v / Pi) 5 ]} 2 - Rotations of R x ,R y ,R z , are car¬ 
ried out for a series of rotation angles as shown in Fig. 
0a). For Raman R x ,R y rotations, the prepared ini¬ 
tial spin state is |s^.). And for Larmor R z rotation, the 
initial spin state is projected by choosing the state of 
idler photon to be | H). In Fig. 0c ~ e), Stokes pa¬ 
rameters are presented for each rotation, and are com¬ 
pared with theoretical curves which fit well. The state 
fidelities of R x ,R y ,R z operations averaged for all rota¬ 
tion angles are 98.6(0.6)%, 98.8(0.5)% and 99.1(0.3)% re¬ 
spectively. Note here the expected state is generated 
by applying ideal operations on the measured density 
matrix of initial states, which takes account of imper¬ 
fections of state preparations and diagnose the rotation 
operations exclusively. To show how well the quantum 
nature is preserved after operations, non-classical corre¬ 
lations between idler and signal photons are shown in 
Fi g- 0g) and(h), where the visibilities are both well 
above 6. Moreover an arbitrary rotation can be decom¬ 
posed as Rn{4>) = R z (ot)R y (/3)R z ( [23] by combin¬ 
ing R y and R z operations. As an example, we choose to 
realize Rn{4>) = exp (-in ■ &<j>) with n = (^=, ^=), 

which is illustrated in Fig. 0b). To make the rotation 
plane meet with zero point, we prepare the initial spin 
state on |i/> s ) = cos|7r|sj.) + sin |7r|s-|-), with a state fi¬ 
delity of 96.6(1.7)%. Then a series of rotation operations 
are performed on the initial state, with roation angles 
4> = 0, y^7t, ..., pj7r. For each rotation, it’s decomposed 
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FIG. 2: Arbitrary rotations of the Spinwave qubit, a-b, II- 
lustratons of state rotations in the Bloch sphere, a, rota¬ 
tions along x-, y-, z-axis. b, special rotation along axis n = 
{^j, 773’ 773}’ an< ^ the rotation plane crosses zero point, c-f, 
Expreimental results of rotations, dotted points are measured 
values of each axis, and compared with theoretical curves, c, 
R x rotations realized by Raman operations with tpR = 0. d, 
R y rotations by Raman operations with ipn = —7r/2. e, R z 
rotations by Larmor operations, f, Ra rotations by combina¬ 
tions of R z and R x rotations, g-h, non-classical correlations 
between idler and signal photons after operations for different 
states of signal photons, g, the initial spin state is on |s^), 
Raman rotations with 1 pn = 0 are applied, h, the initial spin 
state is approximately on |sa) with the idler photon projected 
along | H) state. 


likelihood method (251 . In our experiment quantum pro¬ 
cess tomography is measured for Pauli operations ( a x , a y 
and cr z ) and the Hadamard gate, with the results shown 
in TabpH The process fidelity averaged for these four op¬ 
erations is calculated to be 94.7(7)%. With the process 
fidelities, one can estimate the average state fidelity F ave 
of a process, which is defined as the average of state fideli¬ 
ties between input and output states, where input states 
are randomly selected from the Bloch sphere. It’s shown 
that, the average fidelity of a process, can be simplified to 
evaluate the mean state fidelities of six pure input states 
locating at the cardinal points [25]. The average fidelity 
is related with the process fidelity by, F ave = dFp ^l +1 ■ 
where d is the dimensionality of the system and d = 2 
in our case [27 ■ Both measured average fidelity F ave and 
average fidelity F*% e derived from F proc are also listed in 
Tabjn] We note that F ave is slightly higher than F^ e , 
which is mainly due to the fact that the prepared ini¬ 
tial states are not pure and thus a different definition 
of F(pi,p 2 ) is adopted. Besides, we notice that, the 
constraint of completely positive mapping on the fitted 
process matrixes, also reduces the calculated process fi¬ 
delities, thus further lowering down the F* pe values in 
TabHJ 


T raget 

Operation 

F pmc (°/° ) 

FaveW 

Fdve (%) 

Ox 

<Pr~~ f • t R = 2.6ps 

95.3(1.1) 

97.6(1.6) 

96.9 

Gy 

tp R = 0, t R = 2.6 ps 

93.5(1.8) 

97.7(1.5) 

95.7 


t L = 2.1 ps 

95.3(1.1) 

95.3(3.1) 

96.9 

H 

t L = 2.1 ps 

<PR = -^’t R = 1-3/m 

94.8(1.8) 

98.6(1.3) 

96.6 


TABLE II: Process Fidelity of Pauli and Hadamard Opera¬ 
tions. Pauli operations and Hadamard gate are realized by 
choosing appropriate Raman and Larmor parameters. The 
measured process matrices by maximum likelihood method 
are compared with matrices of ideal operations, to determine 
the process fidelities. Average state fidelities of each operation 
from measured results are also shown here. 


into values of a, /?, 7 , with the global phase term 6 ne¬ 
glected. The result is shown in Fig. if) , with an average 
state fidelity of 98.3(0.8)%. 

In order to give more complete characterization for our 
single-qubit gate, we make use of quantum process to¬ 
mography [24] . An arbitrary single-qubit operation on an 
input state pi n can be fully described by a process matrix 
X, which is defined as p out = Y%,j=o XhjViPin&j, where 
at are Pauli matrices with cro = 7, <j\ = 07 , a 2 = <J y , 03 = 
o z . The distance between paractical and ideal operations 
can be characterized by process fidelity, F proc =tr[xiX'2]- 
In order to measure the process matrix for a single-qubit 
operation, we select six different spinwave states, which 
are listed in Tab|l] as input states and carry out quantum 
state tomography for each output state. In order to get 
a physical process matrix we make use of the maximal 


To summarize, by making use of stimulated Raman 
transition and controlled Larmor procession, we have im¬ 
plemented single-qubit operations for single spinwaves in 
an atomic ensemble quantum memory for the first time. 
We have made use of quantum state tomography to mea¬ 
sure the target states for arbitrary rotations. Average 
state fidelity for R x , R y and R z rotations is measured to 
be 98.8(3)%. We have also adopted the method of quan¬ 
tum process tomography to characterize the Pauli oper¬ 
ations and the Hadamard gate. Average process fidelity 
is measured to be 94.7(7)%. By making use of Raman 
beams with better intensity homogeneity and actively 
controlling the ambient magnetic field, even higher fideli¬ 
ties can be got. Our work enriches the experimental tool¬ 
box of harnessing atomic ensembles for high-performance 
quantum memories, thus could possibly have lots of ap¬ 
plications in future scalable quantum networks [28) . 
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